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Abstract 

Atmospheric parameters, such as pressure (P), temperature (T) and density (p oc 
P/T), affect the development of extensive air showers initiated by energetic cosmic 
rays. We have studied the impact of atmospheric variations on extensive air showers 
by means of the surface detector of the Pierre Auger Observatory. The rate of 
events shows a ~ 10% seasonal modulation and ~ 2% diurnal one. Wc find that 
the observed behaviour is explained by a model including the effects associated 
with the variations of P and p. The former affects the longitudinal development 
of air showers while the latter influences the Moliere radius and hence the lateral 
distribution of the shower particles. The model is validated with full simulations of 
extensive air showers using atmospheric profiles measured at the site of the Pierre 
Auger Observatory. 

Key words: extensive air showers, UHECR, atmosphere, weather 
PACS: 96.50.sd, 96.50.sb, 96.50.sf 



1 Introduction 



High-energy cosmic rays (CRs) are measured by recording the extensive air 
showers (EAS) of secondary particles they produce in the atmosphere. As 
the atmosphere is the medium in which the shower evolves, its state affects 
the lateral and longitudinal development of the shower. Pressure (P) and air 
density (p) are the properties of the atmosphere that mostly affect the EAS. 
An increase (or decrease) of the ground P corresponds to an increased (or 
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decreased) amount of matter traversed by the shower particles; this affects 
the stage of the longitudinal development of the shower when it reaches the 
ground. A decrease (or increase) of p increases (or decreases) the Moliere 
radius and thus broadens (or narrows) the lateral extent of the EAS. 



The properties of the primary CR, e.g., energy, mass and arrival direction, 
have to be inferred from EAS, which can be sampled by an array of detectors 
at ground level. Therefore the study and understanding of the effects of atmo- 
spheric variations on EAS in general, and on a specific detector in particular, 
is very important for the comprehension of the detector performances and for 
the correct interpretation of EAS measurements. 



We have studied the atmospheric effects on EAS by means of the surface 
detector (SD) of the Pierre Auger Observatory, located in Malargiie, Argentina 
(35.2° S, 69.5°W) at 1400 m a.s.l. [1]. The Pierre Auger Observatory is designed 
to study CRs from ~ 10^^ eV up to the highest energies. The SD consists 
of 1600 water- Cherenkov detectors to detect the photons and the charged 
particles of the showers. It is laid out over 3000 km^ on a triangular grid of 
1.5 km spacing [2] and is overlooked by four fluorescence detectors (FD) [3]. 
The SD trigger condition, based on a 3-station coincidence [1] , makes the array 
fully efficient above about 3 x 10^^ eV. For each event, the signals in the stations 
are fitted to find the signal at 1000 m from the shower core, S'(IOOO), which 
is used to estimate the primary energy [5]. The atmosphere is continuously 
monitored by different meteorological stations located at the central part of 
the array and at each FD site. In addition, balloon-borne sensors are launched 
at regular intervals to measure the atmospheric temperature T{h), pressure 
P{h) and humidity u{h) as a function of the altitude h above the detector jlij. 



In section [21 we develop a model of the expected atmospheric effects on 
5'(1000). The modulation is described by means of three coefficients that de- 
pend on the EAS zenith angle {$). They are related to variations of P and 
p, measured at ground level, on slower (daily-averaged) and faster (within 
a day) time scales. The dependence of 5(1000) on P and p implies a mod- 
ulation of the counting rate of events. In section [3l we study the behaviour 
of the recorded rate of events as a function of P and p. On the base of the 
model defined previously, we derive the P and p coefficients. In section |H we 
perform full simulations of EAS developing in various realistic atmospheres 
(based on measurements from balloon soundings above the site of the Pierre 
Auger Observatory) in order to compare, in section [5], the results from data 
and simulations with the predictions of the model. We conclude in section El 
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2 Model of atmospheric effects for the surface detector of the 
Auger Observatory 

2.1 Atmospheric effects on the measured signal 

The water- Cherenkov detectors are sensitive to both the electromagnetic com- 
ponent and the muonic component of the EAS, which are influenced to a dif- 
ferent extent by atmospheric effects, namely by variations of P and p. These 
in turn influence the signal measured in the detectors: for the Auger Obser- 
vatory, we are in particular interested in the effects on the signal at 1000 m 
from the core, 5(1000). 

The continuous measurement of atmospheric P and p is available only at 
ground level. We will show that the variation of 5(1000) can be fully described 
in terms of variation of air pressure and air density measured at the altitude 
of the Observatory site. If not otherwise stated, P and p refer to the values at 
ground level. 

In the following, we flrst describe separately the effects on 5(1000) due to P, 
section 12.1.11 and p, section 12.1.21 and then in section 12.1.31 we provide the 
full parameterisation of its variations as a function of changes in P and p. 

2.1.1 Effect of air pressure variations on the SD signal 

From the point of view of P (which measures the vertical air column density 
above ground), an increase (decrease) corresponds to an increased (decreased) 
matter overburden. This implies that the shower is older (younger), i.e. in a 
more (less) advanced stage when it reaches the ground level. 

The longitudinal proflle of the electromagnetic component of the EAS is ex- 
ponentially attenuated beyond the shower maximum and can be described by 
a Gaisser-Hillas proflle [7] (see Fig. [1]). We are interested in the value of the 
electromagnetic signal measured at 1000 m from the core, referred hereafter as 
Sem- The longitudinal development of the shower far from the core is delayed 
with respect to the one at the core, and can be parameterised as 

5e^(^,X) cxX^— /^exp[(X„,,-X)/A], 

where E is the primary energy, X the slant depth, X^ax = ^max + ^ the 
average maximum of the shower at 1000 m from the core with X^ax being 
the shower maximu A ~ 150 g cm ^ is the typical increase of the shower 
maximum at 1000 m from the core [8] and A ~ 100 g cm~^ is the effective 

^ Xmax — 750 g cm^^ for 10^^ eV showers according to the elongation rate mea- 
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Fig. 1. Average longitudinal profile of three hundred proton- initiated showers with 
E = 10^9 eV, and zenith angle 6* = 60°, simulated with CORSIKA-QGSJETII (open 
blue circles represent the electromagnetic component, red bullets the muonic one). 
The black continuous line is a fit of the electromagnetic profile with a Gaisser-Hillas 
function. 



attenuation length after the maximum [9]. Therefore, a change in P affects 
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where g dX = dP sec 9 is used, with g the acceleration of gravity, and 9 the 
shower zenith angle. Due to the fiat longitudinal development of the muons 
(see Fig. [1]), no significant pressure dependence is expected for the muonic 
component. 



2.1.2 Effect of air density variations on the SD signal 

Regarding p, this affects the Moliere radius tm 

_ EsXq ^ 91m 
~ p ~ p/{kg m-^) 

where Eg = rrieC^ ^ At: / a ~ 21 MeV is the energy constant characterising the 
energy loss due to multiple Coulomb scattering, Ec — 86 MeV is the critical 
energy in air and Xq ~ 37.1 g cm~^ is the radiation length in air. A variation 
in tm affects the lateral distribution of the electromagnetic component of 
the EAS, which can be approximately described with a Nishimura-Kamata- 
Greisen (NKG) profile [TTP^ . At a large distance r from the core, it behaves as 
Sem{r) oc Nem{r) oc T^} {v / T m)'"^ , whcrc T] ~ 6.5-2s and s = 3X/{X + 2X^ax) 

surement with the FD at the Pierre Auger Observatory |T0] 
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is the age of the shower. Hence, a change in p affects 5", 
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In fact, the relevant value of tm is the one corresponding to the air density 
p* two radiation lengths above ground [12j in the direction of the incoming 
shower. This corresponds to ~ 700 mcos^ above the site of the Pierre Auger 
Observatory. On time scales of one day or more, the temperature gradient 
(dT/d/i) in the lowest layers of the atmosphere (the planetary boundary layer, 
which extends up to about 1 km above ground level) can be described by 
an average value of ~ — 5.5°C km~^ at the site of the Auger Observatory. 
Therefore the variation of p* on temporal scales of one day essentially follows 
that of p. An additional effect is related to the diurnal variations of dT/d/i, 
because during the day the surface of the Earth is heated by solar radiation, 
producing a steeper dT/dh in the boundary layer. On the other hand, during 
the night the surface is cooled by the emission of long wavelength radiation: 
dT/dh becomes smaller and even T inversions can be observed before sunrise. 
As a result, the amplitude of the diurnal variation in T (and p) is smaller 
at two radiation lengths above ground than at ground level. It is then useful 
to separate the daily modulation from the longer term one introducing the 
average daily density pd and the instantaneous departure from it, p — pd- 
Therefore, the dependence of Sem on p can be modeled by 



where po = 1-06 kg m is chosen as the reference value of p and is the average 
value measured at the site of the Pierre Auger Observatory over more than 
three years (1 Jan 2005 - 31 Aug 2008). 

Concerning the muonic component of the signal at 1000 m from the core, S^, 
its dependence on p can be parameterised as 



The p dependence is written in terms of pd — po only because, as the muons are 
produced high in the atmosphere, their contribution to signal is not expected 
to depend on the daily modulations taking place in the boundary layer. 

2.1.3 Model of atmospheric effects on S(IOOO) 

The dependence of the total signal at 1000 m from the core, 5'(1000) = S = 
Sem + S^, upon P and p can hence be written as 



em 



5L[i+«r(p^-/^o)+/5r(p-Pd) 



S, = S'^ l + a^^{pd-po) 



S = Sq[1 + ap{P - Pq) + ap{pd - Po) + /3p(p - pd)] 



(3) 
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Fig. 2. Fraction of the total signal induced by the electromagnetic component of 
the shower at ground level at a distance of 1000 m from the shower axis (-Fem) as 
a function of sec 6. A linear dependence of Fem on sec^ (solid and dashed lines) is 
assumed in this work. 



where Pq = 862 hPa is the reference P at the site of the Pierre Auger Ob- 
servatory, 5*0 is the value of the total signal at reference pressure and density 
(P = Po and p = = po), and 

ap = F,^aT «p = i^em«r + (1 - = ^-^r (4) 

where Pgm = Sem/S is the electromagnetic fraction of the signal at 1000 m 
from the core. The values of Pem are obtained by means of proton-initiated 
showers simulated with CORSIKA-QGSJETII (see section H]): they decrease 
approximately linearly with sec 9 for all the simulated primary energies (see 

Fig. ED- 

We will adopt hereafter 

i^e™ = i^:^-0.5(sec^-l) (5) 

where F^^ = Fem{6 = 0) varies between ^ 0.65 at 10"*^^ eV and ~ 0.7 at 
10^9 eV. We note that since the inferred electromagnetic fraction depends on 
the hadronic model adopted and on the CR composition assumed, the actual 
value of Pern may be different. As shown in [U], for iron-induced showers the 
simulated S*^ is 40% higher than in the case of protons, while the SIBYLL 
model [131 predicts a muonic signal 13% lower than QGS JETII for both proton 
and iron primaries. The corresponding variation F^^ at a primary energy of 
10"^^ eV would be ~ —11% for iron with respect to proton, and ~ +4% for 
SIBYLL simulations with respect to QGSJETII. 

Finally, with respect to the coefficients in eq. H) 
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for the pressure coefficient, we have from eq. [T] 

sec 6 



aT ~ -- 
9 



~ X 



A 



where X = X^ sec 6 and X^ ~ 880 g cm ^ is the atmospheric depth at the site 
of the Pierre Auger Observatory. 

(ii) From eq. [2] 

4.5 - 2s 

"p 

where s = 3/(1 + 2cos6' X^ax/X^), with X^ax/Xy ~ 0.85 for 10-*^^ eV pri- 
maries. Pressure effects associated to the change in the slope of the lateral 
distribution function due to the X dependence of s are neghgible. 

(iii) The coefficient /^^"^ should be smaller than a^*" (in absolute value) re- 
flecting the reduction in the amplitude of the p — Pd variations two radiation 
lengths above ground level. The difference should also depend on 9. For in- 
stance, assuming an exponential decrease of the density amplitude with the 
height h 

p{h) - pd{h) = exp [P{0) - PdiO)] 

would lead to 

Z?^'" ~exp(-acos^)a;;'" (6) 

where a parameterises the amplitude of the daily density variation in the 
lower atmosphere and is completely independent of the shower development. 
It characterises the scale height for the decrease of the daily thermal amplitude, 
which becomes 1/e of its ground value at a height (700 m)/a. The value of a 
is expected to be of order unity. 

(iv) The coefficient a'^ is expected to be small, and will be assumed to be 
independent of 9, because of the relatively flat longitudinal development of the 
muons as shown in Fig. [H Its value will be taken to be zero since the air shower 
simulations are consistent with a vanishing a'^ coefficient (see section H]). 



2.2 Atmospheric effects on the event rate 



The dependence of the measured signal on variations of P and p produces also 
a modulation of the rate of recorded events. The trigger probability, Ptr, is a 
well deflned function of the signal [3]. As atmospheric variations correspond 
to signal variations, this implies that the same primary particle (in particular, 
with the same primary energy) will induce different signals depending on P 
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and p. This in turn affects the probabihty for the shower to trigger the SD 
array. 

The effect can be quantified starting from the relation between 5(1000) and the 
energy of the primary cosmic ray. In the case of the Pierre Auger Observatory, 
the primary energy is reconstructed as 

Er oc [5(1000)]-^ , 

where B = 1.08 ± O.Ol(stat) ± O.OA{sys) is derived from the cahbration of the 
SD energy using the FD energy measurement [l^ . Following eq. [3l the primary 
energy Eo{6,P,p) that would have been obtained for the same shower at the 
reference pressure Pq and density po, is related to Er as follows 

Eo = Er[l- ap{P - Po) - a,{pd - po) - Ppip - Pd)f • (7) 

In a zenith angle bin d^, the rate R of events per unit time and unit solid 
angle above a given signal Smm can be written as 

where A is the geometrical aperture and J is the flux of cosmic rays. 

Assuming that the cosmic ray spectrum is a pure power law, i.e. dJ/dEo oc 
Eq'^ , using eq. [71 and neglecting the small energy dependence of the weather 
coefficients, we find that 



ds""^' Is 

oc [1 + B(7 - 1) («p(P - Po) + apipd - Po) + Ppip - Pd))] . 

From the dependence on the atmosphere of the measured CR flux above a 
given signal, we derive the corresponding dependence of the rate of events. If 
Smin is the minimum required signal at 1000 m from the core to trigger the 
array 

— cx [1 + ap{P - Po) + apipd - Po) + &p(p - Pd)\ / dS Ptr{S)S-''^^''-' 

(8) 

with the integral on the right hand side being independent of the weather 
variations. The coefficients ap, Qp and hp are then related to the coefficients 
describing the modulation of the signal by Op^p = P(7 — l)ap,p and hp = 
P(7 - l)/3p. 
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Fig. 3. Top: daily averages of P (left) and p (right). Bottom: diurnal variation of P 
(left) and p (right). The values are averaged over the three years considered (line), 
with the maximum and minimum variations marked by black and white triangles. 
The local time is UTC-3 h (vertical lines mark local midnight and noon). 



3 Atmospheric effects on the experimental rate of events 



To study the modulation of the rate of events, we use data taken by the SD 
from 1 January 2005 to 31 August 2008. All events with 6 < 60° are used, for 
a total of about 960 000 showers with a median energy 6 x 10^'' eV. These are 
selected on the basis of the topology and time compatibility of the triggered 
detectors [Ij. The station with the highest signal must be enclosed within an 
active hexagon, in which all six surrounding detectors were operational at the 
time of the event. 

At the site of the Pierre Auger Observatory, the ground temperature and 
pressure are measured every five minutes. The air density is given by: p = 
(Mm/R) (P/T) where Mm is the molecular mass of air, R the gas constant. 
The daily average density pd is obtained with a smoothing procedure consisting 
in taking, for each time, the average value of p over a 24 h interval centered at 
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Fig. 4. Top: seasonal modulation of the measured (grey) and fitted (black points) 
rate of events. Bottom: diurnal modulation of the measured (grey) and fitted (black 
line) event rate. 

the time of interest. The daily and diurnal variations of the ground P and p are 
shown in Fig. [3] (upper and lower panels respectively). The pressure exhibits 
less than ±2% variation during the period considered, while pd changes up to 
a maximum of ±8% with an additional diurnal variation of density which is 
of ±3% on average with maximum values of 

In the period under study, the number of surface detectors steadily increased 
from about 700 to about 1590. To take this into account, rather than using the 
raw number of triggering events, we compute the rate every hour normalized 
to the sensitive area, which is calculated every second from the total area of 
the active hexagons. The daily and the diurnal rate of events are presented in 
Fig. m (black points), where it is evident that they both follow qualitatively 
the corresponding modulations of pressure and density from Fig. [HJ 

We use the expression given by eq. [H] to fit the measured rate of events. As- 
suming that the number of events rii observed in each hour bin i follows a 
Poisson distribution of average pi, a maximum likelihood fit is performed to 
estimate the coefficients ap, ap and hp. 

The likelihood function is L = H ^e~^\ The expected number of events in 
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bin i is given by 

IJi = RqX AiX d 

where Rq is the average rate we would have observed if the atmospheric pa- 
rameters were always the reference ones, i.e. Rq = ^ , with Ai the sensitive 
area in the i*^ bin and, according to eq. El Ci is 

Ci = [l + ap{Pi - Po) + apiPdi - Po) + bp{pi - pdj]- 
The fitted parameters are: 



ap = (-0.0027 ± 0.0003) hPa"^ 

Up = (-1.99 ± 0.04) kg-^ m^ (9) 
6p = (-0.53 ±0.05) kg~^ m^ 

corresponding to a reduced of 1.06, where = J^ii^i — Pi)"^ / Pi- The result 
of the fit is shown in Fig. |H compared to the daily averaged and the shorter 
term modulations of the measured event rate. 

To check the stability of the coefficients with respect to the energy, the same 
study has been done for the subset of events with a reconstructed energy 
above 10^^ eV, corresponding to ~ 20% of the total statistics. The fitted 
coefficients are consistent within the fit uncertainties. A more detailed study 
of the energy dependence of the coefficients will become feasible in future with 
increased statistics. 



4 Atmospheric effects on simulated air showers 



To complete the study of atmospheric effects, we performed full EAS simula- 
tions in different atmospheric conditions. We simulated proton-initiated show- 
ers using the CORSIKA code [15] with hadronic interaction models QGSJETII 
[ig and Fluka [H]. 

We considered four fixed energies of the primary particle {E = 10^^ eV, 10^^'^ eV, 
10^^ eV and 10^^'^ eV) and seven fixed zenith angles between ^ = 0° and 
6 = 60°. For the air density profiles, we used five parameterisations (shown in 
Fig. E]) of the seasonal average of radio sounding campaigns carried out at the 
site of the Pierre Auger Observatory [6J over a wide range of variation in tem- 
peratur^. The set of simulations consists of 60 showers for each combination 

^ The atmospheric profiles are implemented in the CORSIKA code through the 
dependence of X on h. P, p and T profiles can be derived from: p(h) = —dX/dh 
and P{h) = gX{h). The ground values in Fig. [5] are computed at an observation 



16 




Fig. 5. Left: density profiles used in the simulations. The dashed vertical line corre- 
sponds to the altitude of the Pierre Auger Observatory (1400 m). The corresponding 
values of ground P and T are given in the legend. Right: same density profiles nor- 
malized to an isothermal one = X^/ zq with zq = 8.4 km). 

of atmospheric profile, energy and angle with an optimal statistical thinning 
level of 10-6 [MT9]. 

To compare with model predictions and data, we need to determine for each 
combination {E, 9) the dependence of S'(IOOO) on the variations of P and p. 
The signal can be estimated through simplified assumptions about the energy 
deposited by particles on the basis of their kinetic energy Ek. 

(i) e-e+ deposit Ek — eth, where eth = 260 keV is the energy threshold for 
Cherenkov emission in water. 

(ii) photons deposit E^ — 2me — 2eth- 

(iii) muons deposit 240 MeV corresponding to the average energy released by 
a vertical muon crossing a 1.2 m high water-Cherenkov tank. 

The contribution of each particle is multiplied by the weight assigned by the 
thinning algorithm. We obtain the Cherenkov signal per unit area perpen- 
dicular to the shower plane Csp{r). For the muons, the Cherenkov signal is 
proportional to the track length in the station so that: = C^, whereas for 
the electromagnetic component: C^"^ = cos 6' C|™. 

The left panel of Fig. E] shows the lateral distribution C(r) = C^"^{r) + C^{r), 
which is proportional to S'(IOOO), for four atmospheres (relative to the Spring 
one) in the case of E = 10^^ eV and 9 = 18°. The effect related to the Moliere 
radius can be clearly seen as a broadening of the lateral distribution with 
increasing temperature. 

level h = 1400 m (~ 880 g cm~^), corresponding to the altitude of the Pierre Auger 
Observatory. 
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Fig. 6. Results from the proton shower simulations with E = 10 eV and 9 = 18°. 
Left: lateral distribution of the water Cherenkov signal per unit area perpendicular 
to the shower axis C(IOOO) in four atmospheres normalized to the Spring one. The 
uncertainty is due to shower-to-shower fluctuations. Right: C(IOOO) as a function of 
p for the five atmospheres considered. The dashed and the continuous lines are the 
projections of the fit in the (C(IOOO), p) plane for P = 856 hPa and P = 862 hPa, 
respectively. 

To derive the atmospheric coefficients, we correlate the simulated C(IOOO) 
(taken as the average signal between 950 m and 1050 m) with P and p (see 
eq. [3]). Since we are using seasonal atmospheric profiles, we do not have access 
to the diurnal variation of T and thus we cannot determine the coefficient f3p 
related to the diurnal variation of p. The two coefficients ctp and ap can be 
determined for each fixed energy and angle with a two dimensional fit of the 
C(IOOO), obtained for the five atmospheric profiles, as function of p and P. 
As an example, we show in Fig. [6] (right) the results of the fit for the case 
of E = 10^9 eV and 6 = 18°, projected on the (C(IOOO), p) plane for the 
sake of clarity. Moreover, in the case of simulations we are able to separate 
the electromagnetic and the muonic contribution to the signal and thus to 
determine the atmospheric coefficients for each component (see Fig. [7]). 



5 Comparison among model, data and simulations 

In this section, we compare the atmospheric coefficients derived from data 
with those expected from the model and simulations. We recall that with 
the simulations we cannot access the coefficient Pp, as we use average seasonal 
profiles for the atmosphere, while we can investigate the behaviour of separate 
coefficients for the electromagnetic and muonic components of EAS. On the 
other hand, with experimental data we cannot separate the electromagnetic 
and muonic components, while we can fully investigate the diurnal effects of 
atmospheric changes and compare measurements and expectations for all of 
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Fig. 7. Top: atmospheric coefficients {ap on the left and Op on the right) for the 
electromagnetic component as a function of sec 0. The differently coloured markers 
correspond to the four simulated energies and the lines represent the model for three 
different values oi Xmax- Bottom: ap (left) and Op (right) for the muonic component 
as a function of sec 9. 
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the three coefficients. 



The comparison between atmospheric coefficients for the electromagnetic and 
muonic components of EAS from simulations and model is shown in Fig. [TJ 
as a function of sec 9. With respect to the electromagnetic part, the model 
predictions for both the P and p coefficients, and their dependence on the 
shower zenith angle, are reasonable at all energies. Concerning the muonic 
component of the signal and its dependence on P, ap is compatible with 
zero at all energies, as expected from the flat longitudinal development of the 
number of muons. For the dependence on p, the model is not predictive but 
from the simulations we get a value of compatible with zero. This justifies 
the adoption in the model of vanishing coefficients for the muonic component. 

The comparison of the global coefficients as a function of sec 9 is done for ap, 
Op and I3p in Figs. [8] and [H In the case of the data, the dependence on 9 has 
been studied by dividing the data set in subsets corresponding to five bins of 
equal width in sec 9. For each subset the same ffiting procedure as illustrated 
in section [3] is used. The signal coefficients are then derived by dividing the 
rate coefficients by 5(7 — 1) (see the end of section 12. 2p . Since the bulk of 
the triggering events have an energy < 10^^ eV, we used 7 = 3.30 ± 0.06, as 
measured with the Auger Observatory below 10^^-^^ eV |20] . 

The comparison among data, simulations and model is shown for the pres- 
sure coefficient ap and the daily component of the density coefficient ap in 
Fig. [H] (top and bottom respectively). In the model, we use the value of X^ax 
measured by the Auger Observatory at the median energy of the triggering 
events [10], and a Fg^, corresponding at the same energy, obtained under the 
assumption that F^^ scales linearly with the logarithm of the primary energy. 
The reduced for the data-model comparison is 3.3 for ap and 11.0 for 
ap. For the instantaneous density coefficient /3p, the comparison between data 
and model is shown in Fig. O The data-model comparison gives in this case a 
reduced of 0.6. 



6 Conclusions 

We have studied the effect of atmospheric variations (in P, T and p) on exten- 
sive air showers using about 960 000 events collected by the surface detector 
of the Pierre Auger Observatory from 1 January 2005 to 31 August 2008. We 
observe a significant modulation of the rate of events with the atmospheric 
variables, both on a seasonal scale (~ 10%) and on a shorter time scale (~ 
2% on average during a day). This modulation can be explained as due to the 
impact of the density and pressure changes on the shower development, which 
affects the energy estimator S'(IOOO), the size of the shower signal 1000 m from 
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Fig. 8. Comparison of the P coefficients (top) and of the daily density coefficients 
(bottom) as a function of sec 6 obtained from data (grey shaded rectangle), simu- 
lations (bullets) and model (continuous line). 



the shower axis. This affects the trigger probability and the rate of events above 
a fixed energy. 

The dominant effect is due to the change with the air density of the Moliere 
radius near ground. It induces a variation of the rate of events with associated 
correlation coefficients of ( — 1.99 ± 0.04) kg~^m^ and (—0.53 ± 0.05) kg~^m^ 
on long and short time scales, respectively. 

The second effect is due to the pressure changes, which affect, through the 
variation of the amount of matter traversed, the stage of development of the 
showers when they reach ground. The impact of the pressure variation on the 
rate amounts to (-2.7 ± 0.3) x 10"^ hPa'^ 



21 



" 0.3 

"e 

"a 0.2 
0.1 


-0.1 
-0.2 
-0.3 
-0.4 
-0.5 



Fit of the trigger rate 2005-2008 


^ Modei (F^_^=0.64, X^^^ = 680 g cm"^ 
























1.2 



1.4 



1.6 



1.8 



Fig. 9. Comparison of Pp from data with model. A fit to the data points is performed 
to get the value of the parameter a = 1.7 ± 0.1 (see eq. E]). 



Comparing the coefficients obtained from data, shower simulations in different 
atmospheric profiles and expectations from the model built, a good agreement 
is obtained, not only for the overall size of the effect but also for the zenith 
angle dependence. 

Taking into account the atmospheric effects will allow to reduce the system- 
atics in the energy reconstruction. Furthermore, it will be possible to cor- 
rect for the seasonal modulation, which can affect the search for large scale 
anisotropics. 
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